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Abstract 

We discuss a new approach to three-dimensional electrical impedance 
imaging based on a reduction of the information to be demanded from 
a reconstruction algorithm. Images are obtained from a single mea- 
surement by suitably simplifying the geometry of the measuring cham- 
ber and by restricting the nature of the object to be imaged and the 
information required from the image. In particular we seek to es- 
tablish the existence or non-existence of a single object (or a small 
number of objects) in a homogeneous background and the location 
of the former in the (x, y)-plane defined by the measuring electrodes. 
Given in addition the conductivity of the object rough estimates of its 
position along the z-axis may be obtained. The approach may have 
practical applications. 



1 Introduction 

The aim of electrical impedance tomography (EIT) is to reconstruct the con- 
ductivity distribution o"(x) in the interior of an object Q <Z from electrical 
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measurements on the boundary dQ . For this purpose a number of different 
current distributions are apphed to the surface of the object via electrodes 
and the resuhing potentials on the surface are recorded. Applications can be 
envisaged both in medicine and industry 

Conservation of the current j(x) and Maxwell's equations in the quasi- 
static limit lead to the following differential equation for the potential $(x): 



In the following we take as the object a rectangular box and investigate 
whether statements on the conductivity distribution can be made if the sur- 
face potential can only be measured on one side of the box. Such a model 
relates to typical situation in geological and medical imaging. 

The general inverse conductivity problem for the box requires current- 
and potential-measurements for a large number (in principle infinite) of ap- 
plied current configurations on the surface of the box. For the reconstruction 
of the conductivity distribution in this and related problems the boundary 
conditions must be known precisely and all calculations of potentials be per- 
formed with high accuracy. All these conditions are difficult to be achieved 
in practice, which explains the comparative lack of success of the impedance 
method in medical applications. In many cases, specifically breast cancer 
screening, it is actually not absolutely necessary to have a complete image of 
the region. If we restrict the reconstruction to a shadow on a plane and re- 
quire only rough information on size and location of the cancerous region, the 
reconstruction can be done analytically using a single measurement. This 
problem has also been discussed from different points of view 0. 

2 Description of the problem 

We are interested in the conductivity distribution (t{x) inside a rectangular 
box with sides a, b, c, as pictured in figure (|Tp. 

The region of interest ft is therefore of the form 





Q = {{x,y,z)\ 0<x<a, 0<y<b, 0<z<c}. 



The boundary is made up of six rectangles 



dQ = dQ 



'1=0 



UdQ. 



'x=a 



u9a 



'y=o 



UdQ. 



•y=b 



•Z=C ) 
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Figure 1: The geometry of the imaging device 

where, for instance, dQx=o means 

dnx=o = {(x, z)\x = 0, < y < b, < z < c} . 

Similar definitions hold for the other rectangular regions. 

The following discussion assumes that a fixed external current enters on 
one of the side surfaces and leaves on the opposite surface. The current is 
taken to be constant on the two surfaces, i.e. 

(x, y, z) e dnx=o (2) 
(x, y, z) e dVtx=a (3) 
otherwise, (4) 

where ^ denotes the normal derivative. For simplicity we set /= 1, which 
can always be achieved by a suitable choice of units for the conductivity and 
the potential. 

We further assume that conditions are such that the resulting potential 
can only be measured on the plane dVtz=Q (see fig.(|l])). 

Given cr(x) the resulting potential 0(x) can be obtained by solving the 
differential equation (P with the Neumann boundary condition Eq.(H,||,^. 

The aim is to obtain an image of cr(x) from the measurement potential 
on the boundary dVLz=Q. If the conductivity does not differ much from a 
constant distribution ctq, we can write 

O-(x) = (To + ^cr(x). (5) 



d(j) 



dn 
c 

a 



dn 



-I; 
= /; 



^7r = 0' 

on 
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Without loss of generality we can set ctq = 1. For a = ao the solution of the 
boundary value problem is obviously 

(j)o{x,y,z) = X + const. (6) 

In the following section we try to answer the question to what extend 
5cr(x) can be reconstructed by measuring the potential only on the lower 
surface of the box, i.e. on the boundary surface SQ^^o- 



3 Reconstruction 

As the potential distribution is only defined up to a constant, it is conve- 
nient to require that the average of the potential distribution vanishes on the 
boundary surface dfl 



= (7) 

on 

If we assume that the current on the surface dQ is square integrable, it 
is in the space 



ClidQ) := !^feC\dQ), = o}. (8) 

Any change 6a of a homogeneous conductivity distribution ctq produces 
a corresponding change 6(f) in the potential distribution 0o- Then, for any 
function g G Cl{dQ), it can be shown (see appendix), that in linear approx- 
imation 

< H,9 >c'i(dn)= / S(pg = - V(pQ-V(pg6a , (9) 
Jan Jn 

where 0^ represents the solution of the potential problem for constant con- 
ductivity o"o = 1 and external current distribution g [Q. We have checked 
in model calculations that the linearization yields good qualitative images 
even for objects with large conductivity. This refers only to the geometri- 
cal appearance of the objects and not to the actual numerical value of the 
reconstructed conductivity. A small spherical metallic object of infinite con- 
ductivity in a homogeneous background of conductivity 1 unit, for example, 
will be imaged as an object of conductivity 3 units. We will therefore not 
attempt to determine numerical values of the conductivity of the hidden ob- 
jects but only existence and location of such objects. This restriction of the 
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scope of the reconstruction will still yield useful results in applications such 
as mammography. 

For the given experimental set-up we measure a change in potential 6(j)exp, 
which we normalize so that it is in Cl{dQz=o), which is defined in analogy 
to Eq.dl). We consider a base {un}, which is complete and orthonormal in 

It turns out to be useful to introduce in addition a set of functions m„ G 
Cl{dQ) (not complete), which are defined on the full surface dQ of the box, 

-(^) - { ""'1 ; • (1°) 

Then, by Eq.(|^), the moments < Scpexp, Un >cl(dn^^o) satisfy in linear approx- 
imation 



- [ V0O- V0^„ 6a. (11) 
Jn 

We introduce a linear operator A acting on the change in conductivity 
Sa through 

^<5a:=-^ V0O- V0^„<5a^ Un. (12) 

Using Scpexp = J2n ^ ^(pexp, Un >c'i{dn^=o) relation between 6a and the 

associated change in potential 6(f)exp reads 



A5a = S(P,.,p. (13) 
A natural choice for the base {u„} associated to the upper surface is 

(z,j)^(0,0), (14) 
(15) 

where the index n is replaced by two indices The set of functions 
Un — ^ Uij referring to the whole surface is then defined in accordance with 
Eq.(0). To make use of Eq.(|D to calculate 6a we need the potential (pa. , re- 
sulting from an external current distribution Uij G Cl{dQ) and conductivity 
(To = 1. It is a simple exercise to show that 

0..,, (X, y, z) = ^^^^_^,s.^^^ cos — cos — {e + e } , 

(16) 



Ui,j{x,y) 



l-KX jTxy 
Lij cos cos 



a 

2/v/a6 
v/27p) 



h 
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with the abbreviation 



= n^{z/ay + {j/by. 

If we define 



aij = ||V0o ■ 



i TT 



adij{l - e ^-^'jc^ ^^^^^^ 



1/2 

(17) 



-V0O ■ 



||V0o ■ V0n^ j£2(f^) 

C. . sin — cos ^ {e-^^-^ + e^'-^-^^'-n 

-1/2 



then Eq.([T^) can be written in the form 



AScr = ^ (Tij < S(r,Vij >c2(n) Uij- (19) 

i=l,j=0 

This is our main resuh. It is obvious from Eq.(|TU]) that the set {vij} is a 
complete orthogonal system in N{A)-^. We have thus explicitly constructed 
the singular system {vij,Uij;aij} of the operator A and its generalized in- 
verse can be written down explicitly. The generalized or least square solution 
of Eq. (|19D is then simply given by 

oo 

6a = ^ aji'j < 6(pexp,Uij > Vij. (20) 

i=l,j=0 

This generalized solution is still not continuous in the data and must be 
regularized in a suitable manner. We employed for convenience mainly the 
method of truncating the singular values or truncating the indices i,j in 
Eq.(^). The latter procedure turns out to produce better images. The 
cut off values of the indices is determined by a version of the discrepancy 
principle, i.e. by requiring that the resolution implied by the Fourier series 
Eq.|l^ should not exceed the distance between the electrodes which measure 
the potential on the surface dflz=o- We assume the latter constitutes the 
main source of the experimental error. 
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It should be pointed out that we only reconstruct a three dimensional 
picture which is a projection on the set of functions {vij}. It is sufficient to 
view the image at z = 0, because the images for z follow uniquely from 
the z = one and contain no additional information. 

We effectively see a two-dimensional image, which represents a kind of 
shadow of the object. For many purposes (such as in cancer screening), when 
one is only interested in the presence or absence of an object, this is sufficient 
information. As discussed above, the actual value of the object's conductivity 
cannot be reconstructed quantitatively. 

In Fig.(|^) and (0) we present images obtained from synthetic data which 
were calculated for 10 x 10 grid points on the plane dQz=o ■ We also show 
how the image deteriorates when errors are assigned to the data. 



a b 




024 68 10 024 68 10 



Figure 2: Images of an spherical object obtained from exact and error affected 
data 

Figure @ shows images of a spherical metallic object of diameter d = 1 
at a distance z = 2 (all in units of the grid spacing) from the surface of 
measurement, (a) with exact data, and (b) with data 6(j)exp corrupted with 
a 20% random uniform multiplicative error. It is amazing that even with 
errors of such a magnitude a reasonable image is produced. 
Figure (||) shows images of two spherical objects obtained from exact data. 
Case (a) shows the image for two spheres of diameter 1.5 and 1 respectively 
both at a distance z = 2 from the surface and case (b) shows the image for 
two spheres of equal diameter at distance z = 1.5 and z = 2 respectively. 
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a 



b 
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Figure 3: Images of two spherical objects 

Qualitatively the image gets larger and flatter as the object is moved 
away from the measuring plate. In addition the image gets brighter (but 
not larger!) as the object gets larger. The same effect is observed when 
the conductivity is increased. It is not possible to distinguish volume- from 
conductivity effects. This is true as long as the objects are small or not 
to close to the surface, as they effectively behave as dipoles (see section 4). 
Given additional information, e.g. that the object's conductivity is constant 
and of a given magnitude, it may be possible to quantify this observation 
and obtain a full three dimensional image of the object. This is exemplified 
in the next section. 

4 Spherical Object 

In the following we consider a single spherical object K of conductivity n 
and radius a immersed in the box Q filled with a liquid of conductivity 1. 

Let uk be the normal to the surface of K and uq the normal on 
The boundary value problem for a current distribution / e Cl{dQ) can be 
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defined as follows, 





A(j){x) 


= ,x 


en\dK, 


(21) 




d(f) 
dnn 


= /(^) 


, X G dil, 


(22) 


lim (dlx + hnx] 


1 — (f){x — hriK)) 


= ,x 


e dK, 


(23) 


/ d4>{x + hrix) 
V driK 


d(j){x — hrix) \ 
driK ) 


= ,x 


e dK, 


(24) 




I (pds 
Jan 


= 0. 




(25) 



lim 



Equation (^) guarantees the continuity of the potential while (^4D describes 
current conservation. These two equations determine the boundary condi- 
tions on the surface of the sphere. The other three equations represent the 
well-known boundary value problem of the Laplace equation. The Neumann 
boundary condition ( P^ is given in (0,^,^). For k = 1 one obtains (po, the 
solution of (^. 

For the case of a small sphere of constant conductivity and not too close 
to the surface Q, the change in potential 6(j) is given by the dipole term, 

V0O • jx-xp) 

50(x) = — — ; (26) 

\x — Xol 

where Xq is the coordinate of the centre of the sphere, and 



This result can be derived by noting that our boundary value problem is 
equivalent to that of a dielectric sphere in a uniform electric field. 

The variation of a with k shows quite clearly the limited sensitivity of 
EIT to changes of conductivity. In practical reconstructions k = 10 can 
hardly be distinguished from k = oo. 

The potential (p = (po + S4> still does not satisfy the Neumann boundary 
condition ( p2D on dQ. This problem can, in principle, be solved by an infinite 
number of image dipoles. We checked that the series converges rapidly. For 
the case a,b ^ oo , i.e. the case of two infinite plates,the sum takes on the 
simple form: 
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OO f 

5(j)[x) = -ar^ ^ < ■ 



V0O ■ ix-Xo,y- yo, z - {2nc - Zq))' 



„=o V ,(3; - 3^0)^ + {.y- yoY + {z- {2nc- zo)y\^ 
V(f)o- {x -Xo,y -yo,z - {-2nc + Zq))'^ \ 
^ \(x-xoy + {y-yoy + {z-{-2nc + zo)y\'j ^ ^ 

The knowledge of 6(f) allows to calculate the generalized inverse according 
to Eg ■ (|20|) ■ For rough estimates and when the sphere is not to close to the 
surface, the image dipoles can be neglected. As the conductivity of the sphere 
is assumed to be known, Eq.(^) can be used to obtain an estimate of the 
position and the volume of the sphere. As a test we measure the synthetic 
data on the surface dQz=Q by 10 x 10 electrodes. Given the rough knowledge 
of the coordinates (xq, yo) of the centre of the sphere, we fit 6(f) of Eq. (|26|) plus 
the background potential 0o to data taken on neighboring electrodes. As a 
typical example, we find for data afflicted with 10% multiplicatice uniform 
error, the following results: Zq = 2.32 (instead of 2.0) and r = 0.56 (instead 
of 0.5). 

In realistic applications the object to be detected will in general not be 
spherical. Nonetheless one may obtain rough information on size and depth 
of the location of the object by assuming a spherical shape and applying the 
analysis above. 



5 Conclusion 

We have presented in this note an electric impedance imaging system based 
on a specific simple geometry of the device which guarantees a uniform cur- 
rent distribution in the case of constant conductivity. If we further impose 
the condition that only a single object (or possible a small number of ob- 
jects) is to be detected, then we show that an image can be obtained in 
a single measurement of the surface potential. To test of the effectiveness 
of the method, we create synthetic data which can be afflicted with errors. 
The image obtained by inverse problem techniques represents a projection or 
shadow on the surface where the potential is measured. This image is amaz- 
ingly stable against data errors. We also indicate how rough estimates on the 
size and the depth of the object may be obtained. The actual construction 
of such an imaging system is planned. 
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A Appendix 

We will give the sketch of a proof of Eq.(P). Let a potential, denoted by 
m(x), satisfy the EIT differential equation 

V((tVm(x)) = (29) 

for a given conductivity o"(x) and surface current /(x) 



f - a 

On 



(30) 

an 



Let f (x) be an arbitrary solution of the EIT differential equation Eq.[29|. 
Then we can define a functional 



hj[v]:= f-v (31) 
Jon 

and a bilinear form. 

a„[u,v] := / (xVn ■ Vv . (32) 
Jn 

The EIT boundary value problem with Neumann boundary conditions is 
known 0]to be equivalent to the condition 

bf[v] = a(j[u,v] Wv (33) 

We now change a a + 6a, keeping the same Neumann boundary condition 
30| . Then the potential will change as u ^ u + Su and the condition Eq.p^ 
will read 

bf[v] = a^+sa[u + Su, v] (34) 

Or 

bf[v] = a^[u,v] + a^[6u,v] + asa[u,v] + asa[Su,v] (35) 

Neglecting the last term,using Eq.^ and the symmetry of the bi-linear form, 
we obtain the relation 

asa-[u,v] = —ao-[v,6u] Wv (36) 

As f (x) is the solution of Eq.|2^ for some boundary current 5'(x), it must 
satisfy 

bg[w] = a^[v,w] (37) 
for all w(x), in particular for w{x.) = 6u{x.). We finally obtain therefore 

g ■ 5u = —asa[u,v] , (38) 



an 



which is just Eq.(9). 
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